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Abstract 
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The random networks enriched with additional structures as metric and group- 
symmetry in background metric space are investigated. The important quantities 
like he clustering coefficient as well as the mean degree of separation in such net- 
works are effectively computed with help of additional structures. Representative 
models are discussed in details, 
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1 Introduction 

It is well-known for about 40 years that the mean degree of separation of humans 
(the mean number of people needed to bind two randomly chosen persons via a 
chain of acquaintance) is surprisingly small, it takes the value about six. This was 
discovered for the first time by the American psychologists Stanley Milgram pQ. 
Since his days this very interesting phenomenon attracts the interest of scientists 
from various fields of research. As shown in [6], the small world phenomenon 
appears as a generic feature of many natural as well as artificial random graphs. 
To explore these properties of the random networks and to understand them one 
finds often helpful to use methods developed for statistical mechanics [5]. Gen- 
erally, there is no hope to obtain exact results for main quantities describing 
random networks of real interest. However, under some additional condition one 
can build up simplified models obeying some kind of symmetries that allow for 
expressing crucial quantities (their approximations) in closed form. Methods of 
computation can be divided into two classes; one of them uses purely discrete 
point of view and the second one, which is much more effective if applicable, uses 
continuous limit of the network. We shall be oriented in the latest. 
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The importance of the use of a continuum limit in various discrete random struc- 
tures is surely undisputable because the exact discrete mathematics based meth- 
ods have commonly bounded area of use. In current research the continuum limit 
methods are used widely, for the use in direct random graph theory see e.g. [2J, 
an interesting result is achieved by these methods studying some hidden variables 
(variables that are not measured directly like degree of a node) in a kind of bi- 
ological networks [3]. One of the most important areas of research where these 
continuum limits are used is the problem of phase transition in random structures 
(critical phenomena), for the review see e.g. [3] and further references therein. 
The main aim of the present paper is to show how to compute the relevant quan- 
tities like the clustering coefficient and the mean degree of separation of nodes in 
a given random network that is build up on a metric space or a space where we 
can define some oriented distance. The crucial assumption is that the network 
obeys some symmetry. The metric serves us to define the probability that two 
points (nodes, persons) are linked via a kind of acquaintance, and this probability 
is a given function of the distance between the nodes. 



2 Definitions and basic properties of relevant 
quantities 

Let us consider the number N of entities (people) represented as points lying in 
the metric space f2 equipped with the distance d. The acquaintance (defined in a 
suitable way) between two people represented by the points A and B is recorded 
as a link between these points. With the help of the metric considered on Q we 
can construct a probability Q(A, B) that A knows B (and vice- versa, we consider 
always the reflexivity of knowing someone) as a function of the distance between 
A and B: 

Q{A,B) = Q{d{A,B)). (1) 

The simplest quantity characterizing our network that is given by Q is the average 
number of acquaintance for any person A: 

m=j2q(ab), (2) 

B 

where the summation runs over all points of considered network. If one assumes 
that the function Q varies slowly on the scale of elementary distance in the 
network, then the above written sum can be approximated by the integral: 

N= [ Q(x,y)dy (3) 
Jn 

with proper measure dy on Q. (For example, Q can be considered to be an sub- 
manifold in some Euclidean space M n and then it is equipped naturally with both 
metric and volume form.) 
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Let b be the distance between the two considered points. Then the average 
degree of separation can be written as the mean value of the function (mean 
with respect to its argument S) P(S, b) that is the probability that the degree of 
separation of A and B at relative distance b equals S. 

First of all, P(0, b) is given trivially by definition of Q itself: 

P(0,b) = Q(b). (4) 

Other P(S, b)'s are to be computed nontrivially. And we are able to obtain ana- 
lytical results only if some approximations in computing P(S, b) are applicable. 
Let us discuss P(l, b). This means that there is just one person, say C separating 
the two chosen persons A and B. (Let us note that C is necessarily different from 
A as well as from B, otherwise A knows B.) With respect what has just been 
said we obtain the probability in question: 

P(l, b) = J2 Q(A, C)Q(C, B)(l - Q(A, £?))• (5) 

c 

Analogically we comes at the expression for the probability that the separation 
of the two points at fixed distance is two: 

P(2, 6) = Y, E Q( A > Ci)Q(Ci,C 2 )Q(C 2 , B)(l - Q(A, C 2 )) x 
G\ C2 

(l-Q(B,C 1 ))(l-Q(A,B)), (6) 

and we could continue to further separation indices in the same way. However, 
instead of doing this, we are going to describe an approximation we shall use in 
the future considerations. We assume that the probabilities Q(A, B) are small 
enough and therefore the product of two such probabilities (that is of the higher 
order of smallness) can be neglected in the expression ([6]) as well as in others 
P(S,b)'s. In this way one obtains instead of JS]) its reduction: 

P(2, b) « E ^ C M C i> C 2 )Q(C 2 , B). (7) 
Ci c 2 

By means of derivation ([3]) from ([2]) we can write down for P(2, b) also the fol- 
lowing integral formula 

P(2,6)« 

Q(A, x)Q(x, y)Q(y, B)(l - Q(A, y))(l - Q(x, B))(l - Q(A, B))dxdy 

nxn 

« ( Q(A,x)Q(x,y)Q(y,B)dxdy. (8) 
Jnxn 

Obviously, the same can be done for all other P(S, 6)'s with the similar result. 
The special structure of the integral formula ([8|) will allow for interesting expres- 
sion for the essential quantity we shall turn now our attention to. The quantity 
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that is effectively and often used in analysis of random networks is the mean 
clustering coefficient (C), [7J. The quantity (C) is the probability that (any) two 
acquaintances of given person A knows each other. By this definition, (C) can 
be computed modifying equation ([5]) - namely we must identify A and B and 
avoid interchanging the positions of persons C\ and C2 that means to reduce the 
average number of triples by the factor 1/2: 

- / Q(A,x)Q{x,y)Q{y,A)dxdy. 

To obtain (C) we have to divide the above written average number of triples by 
the number of possible triples that is given by the fact that average number of 
acquaintances for any person is J\f, so the normalizing constant we are looking 
for reads: N{J\T-t) ^ hP_^ where we have assumed H 3> 1. Finally, the clustering 
coefficient is given by the following simple equation: 

^ C ) = Tnl Q(A,x)Q(x,y)Q(y,A)dxdy. (9) 

The idea of the following is that the above written quantities can be expressed 
in approximate closed form using additional assumptions on symmetry of the 
metric space in use. We shall demonstrate this in details on a concrete example, 
however, the idea can be extended to wide range of symmetric spaces in the same 
way as generalized Fourier series are introduced. 



3 One dimensional closed model 

Let us consider the standard circle S 1 (one dimensional sphere) as our metric 
space Q. This space can be considered as the orbit of a point under the action of 
the U(l) (or alternatively SO (2)) group. The position of any person is then given 
by its polar coordinate <p ((f) € [— 7r,7r]) only. The probability Q is the function of 
one real variable by means of 

Q(<f>u<h) = Qifa - fa) = - (10) 

The topology of the problem tells us that Q is a period function with period 2tt. 
This implies we can expand the function Q into its Fourier series in trigonometric 
functions in the mentioned interval [— ir, 7r]: 

00 

q{<i>) = y, a ^ lH - 

k=—oo 

Taking into account ()10p (and that Q is real-valued) we have that the Fourier 
coefficients obey the equation: a_fc = afc- As an example we shall compute the 
P(l,b) and P(2,b) functions in our case. With respect to eq. ([S]) and using 
orthogonality property of trigonometric basis we have after some algebra 
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P(l,b) = 2irR(l-Q(b)) ^ a 2 m cos(mb) , (11) 



m=— oo 



where the parameter R is introduced as the radius of the considered circle. This 
parameter will not enter the results because its serves only as a geometrical tool 
to keep the distance between any two neighbours equal to 1 if suitable. This 
simply means that the volume form is Rdcj). Supposing the persons A and B are 
distanced enough (6^1 => Q(b) <C 1) we can write instead of eq. (|lip its 
simpler form: 

oo 

P(l,&)«27rR Yl a^cos(m6). (12) 

m=— oo 

Using the same idea we can find also the probability P(2, b): 
P(2,b) = 

( oo 

{2,fR\l-Q{b))\ £ ale—" -2 £ £ 



oo oo 
m=— oo n=— oo 



oo oo oo 



^ ^ ^ ^ ^ ^ Q / m+n+pQ'rn+'na"mQ'nQ'p& ^ ^ I (^) 
ra=— oo n=— oo p=— oo ) 

This expression is a little bit more complicated, assuming, in the same way as 
previously, b is large enough and neglecting the sums consisting of higher than 
third power of the Fourier coefficient we have the approximative expression for 
the probability P(2,b): 

oo oo 

P(2,b) « (2tt) 2 R 2 Yj a 3 m e- imb = (2tt) 2 R 2 ^ «&cos(ro&)- (14) 

m=— oo m=— oo 

Thus, the important quantity as the clustering coefficient is expressed as follows 
immediately from eq. (|13p putting 6=0: 

4^2^)2 r 00 oo oo 

^m=— oo m=-oo n=-oo 

oo oo oo \ ^2 R 2 °° 

m=— oo n=— oo p=— oo J m=— oo 

Furthermore, one could easily generalize the results (|12p and (I14p to the proba- 
bility P(k,b), the result reads: 

oo 

P(k,b) = (2nR) k a* +1 cos(nb). (16) 

n=— oo 

Now we shall apply our theoretical results to a concrete example of primary prob- 
ability function Q often discussed in literature in this context. 
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Uniform distribution within a fixed radius 



Now, the function Q will represent the uniform distribution of probability 
within a fixed radius (angle), namely: 



Q(<f>) 



p , 

, <££[-$,$] ' 



(17) 



where $ is a parameter lying in the interval [0, ir] and p is the probability of an 
acquaintance, therefore p € [0,1]. The mean number of acquaintances of any 
given person in this model is given by 



TV = 2R / pd(j> = 2Rp$. 
Jo 



In order to obtain the clustering coefficient we are to find the Fourier coeffi- 
cients of the function (|17p in the interval [— 7T,7t]: a n = J- j^_ iT Q{cj))e m< l > dcf) = 
HoQW cos(n^)d0; 



a 



7T 



pdtfi 



p<3? 



7T 



* / psin(n<I>) 

pcos(n<p)a<p — 



o 



7T n 



With respect to fj 1 5|) the clustering coefficient is given by: 



(C) 



7T 



$2 



$ 3 + 2]T 



sin 3 (n$) 



n=l 



(18) 



Special case arises when 3> = tt, i.e. when one considers uniform distribution of 
probability throughout whole circle. In such a situation we have trivial result: 



(C)($ = 0)=p. 



(19) 



The result (|18p is plotted in the graph [TJ The probabilities P(k,b) in this model 
are given by (fT6|) : 



P(k,b) = (2nR) k 



n=l 



psin(n$) 



irn 



k+1 



cos(n6) 



p /AT 

7T V $ 



ra=l 



(20) 



and for the maximally distanced persons: 



P(fc,7r) 



p /A/" 



7T V $ 



n=l 



(21) 
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Figure 1: Left graph shows the plot of the fraction (C)/p for the clustering coefficient 
given by ()18|) . The fraction is nearly a constant (except the value at $ = that is, 
however, trivial) up to an angle a little bit greater than 2 and then starts to rise up 
to 1 as it should with respect to f[T^]) . Right one shows the plot of the functions: 
P(k,7r) = P(k,n)/(7iAf k ) in dependence on the angle The curves are plotted (from 
the right to the left) for k = 1, 2, 4, 6, 10, 20. One can obviously identify critical values 
of the ^-parameter for (C)/p as well as for P(k, ir). 




4 Higher dimensional models 

Now we shall briefly show how the previous consideration can be extended into 
higher dimensional case. For example, people living on Earth are to be modeled 
as, at least, two dimensional network (in fact one should incorporate also socio- 
logical dimension rather than geographical only). 

There is an straightforward possibility how to generalize our one-dimensional 
example. Instead of S 1 we can consider the iC-dimensional tori T , i.e. the 
network with the topology 

T K := S 1 x S 1 x • • • x S 1 . 

* v ' 

X-times 

The edges between any two neighbours have the unit length. Let (</>i, 02, • • • , <Pk) 
be the (periodic) coordinates. The probability that the persons A = (0, . . . , 0) 
and B = (0i,...,0#) knows each other is given by the functions of angular 
distances between the persons in each direction: 

Q(A, B) = Qx^x) • Q 2 (<p 2 ) ■ ■ ■ Qk(<^k), (22) 

where any function Qi plays the role of single function Q from the previous 
section. The prescription (|22p for the probability Q allows for introducing an 
anisotropy in the model, if suitable. This situation in two dimension is shown in 
figure [2 Let us mention that this example shows we do not need to have metric 
space, the metric structure is replaced in this case by the capability to measure 
distances in separate directions. 
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The task to find the probabilities P(k, b±, ... , bx) can be obviously solved in 
the same way as it was done in one dimensional case in previous section. The 
difference stands in only one detail that the single integrals (defining Fourier co- 
efficients) are replaced by multiple (A-tuple) integrals defining the Fourier coeffi- 
cients of functions defined on S 1 x S 1 x • • • x S\ . Of course, a practical difference 

X-times 

can appear, namely it is reasonable to expect the numerics in higher dimensional 
case will be more complicated. Moreover, one can use the same procedure for 
any homogenous space with help of generalized Fourier series. Especially, one can 
perform analogical computation also in the case when (at least in one dimension) 
the background space is unbounded - the only thing is to replace Fourier series 
by Fourier (generalized) transformation. Generalized orthogonality will allow for 
analogical expression of the clustering coefficient also in the case of non-compact 
symmetry group. 

flat torus T 2 = S 1 x S l 




Figure 2: Two dimensional flat torus as two dimensional compact background space. 
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